Second-order equation of state with the full Skyrme interaction: toward new effective 

interactions for beyond mean-field models 
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In a quantum Fermi system the energy per particle calculated at the second order beyond the 
mean-field approximation diverges if a zero-range interaction is employed. We have previously 
analyzed this problem in symmetric nuclear matter by using a simplified nuclear Skyrme interaction, 
and proposed a strategy to treat such a divergence. In the present work, we extend the same 
strategy to the case of the full nuclear Skyrme interaction. Moreover we show that, in spite of the 
strong divergence (~ A 5 , where A is the momentum cutoff) related to the velocity-dependent terms 
of the interaction, the adopted cutoff regularization can be always simultaneously performed for 
both symmetric and nuclear matter with different neutron-to-proton ratio. This paves the way to 
applications to finite nuclei. 



PACS numbers: 21.60. Jz,02.60.Ed,21.30.-x,21.65.Mn 



I. INTRODUCTION 



Mean-field theories are widely employed to analyze a large variety of many-body systems where an independent- 
particle picture can be adopted. In nuclear physics, mean-field approaches generally lead to satisfactory results when 
applied to several bulk properties of atomic nuclei such as masses, radii, or ground-state deformations. In some cases, 
however, in nuclear physics as well as in other domains of many-body physics, mean-field approaches need definite 
improvements. For instance, in the nuclear case, mean-field models do not predict accurately the single-particle 
spectra (namely, energies, spectroscopic factors and fragmentation of the single-particle states, especially close to the 
Fermi energy): the mean- field Hartree-Fock (HF) occupation numbers are strictly equal to or 1 and, therefore, the 
spectroscopic factors of single-particle states cannot be reproduced. To improve the theoretical description, one may 
introduce the coupling between nucleon individual degrees of freedom and collective degrees of freedom by means of 
the particle- vibration coupling approach [IHIjj. This model constitutes an example of beyond mean-field theories. 

A conceptual problem is often encountered when beyond mean-field theories are employed in nuclear physics where 
the most widely used interactions are phenomenological. That is, the parameters of these interactions already take 
into account in an effective way several correlations. Hence, their use in beyond mean-field models (where correlations 
are explicitly introduced) may lead to a double counting. This difficulty could be removed in principle by adjusting 
£S) ' the parameters directly at the beyond mean-field level. 

t-H ' An additional practical problem appears in beyond mean-field theories where the interaction has zero range. In the 
^ . nuclear case, all the terms of the Skyrme interaction, and the density-dependent and spin-orbit terms of the Gogny 

• rH ' interaction have zero range. Due to this, ultraviolet divergences are generated in the diagrams beyond HF, so that 
, the numerical results are dependent on the choice of the energy cutoff. 
H 1 Different procedures to regularize the divergences that appear beyond HF may be adopted. For example, in the 
context of the Hartree-Fock-Bogoliubov or Bogoliubov-de Gennes model (where an ultraviolet divergence is generated if 
a zero-range pairing interaction is used) a currently employed regularization is based on the so-called pseudo-potential 
prescription This procedure has been proposed in nuclear [f| and in atomic @, 0] physics. Another example is 
the cutoff regularization which is adopted in this work, where we change the cutoff and adjust simultaneously the 
parameters of the Skyrme effective interaction to keep the energy per particle at the beyond mean-field approximation 
independent of the cutoff regulator. A third example is the elegant procedure of the dimensional regularization which 
has been developed in the 70s in the context of the electroweak theory [§4lbT| and which is widely employed nowadays in 
perturbative quantum field theories to regularize the divergent integrals. The basic idea of this regularization technique 
consists in replacing the dimension d of the integrals with a continuous variable. The dimensional regularization (at 
variance with the cutoff regularization) preserves translational, gauge and Lorentz symmetries [TTj and this explains 
its wide diffusion and success in the framework of quantum field theories. 

The ultraviolet divergence appearing beyond the mean- field approach has been discussed in a previous work 
where the nature of the divergence has been analyzed in the case of symmetric nuclear matter within a simplified 
to — t% Skyrme model. The used interaction was V{r\,r2) = g{p)5{f\ — r^), where the coupling constant is density 
dependent and equal to g(p) = to + ^p a \ to, ^3 and a are parameters. The second-order correction to the mean- field 
equation of state (EoS), that is, the lowest-order term beyond HF, has been studied analytically in detail. The precise 
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nature of the divergence, which is linear in the momentum cutoff A, can be put in evidence by using the results of 
Ref. [l2j . In such a case, the linear divergence in A can be explicitly shown by taking the asymptotic expansion (for 
A — > oo) of Eq. (8) of Ref. [ijj. This expansion reads 

-ll + 2]og2 , A 2k F # 
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For different values of A, new sets of parameters (to, t 3 and a) have been adjusted to reproduce a reference EoS 
for symmetric matter. In practice, the mean-field EoS of the SkP Skyrme force [13| | has been chosen as a benchmark. 
This choice has been suggested by the fact that, for the SkP paramctrization, only the to and t 3 terms actually appear 
in the EoS for symmetric matter. In the framework of such cutoff rcgularization, the choice of A has been dictated by 
physical criteria. As already mentioned in Ref. fl2j , since in low-energy nuclear physics nucleons are treated as point 
particles, A must be definitely smaller than 2 fm . This value of A will be the maximum cutoff adopted in this work. 

In the present work, the approach of Ref. [l2j is considerably extended. The full Skyrme interaction is taken into 
account and both symmetric and asymmetric nuclear matter (including the limiting case of pure neutron matter) 
are considered. The latter point is particularly important since we show here that the obtained interaction fitted 
for several values of A to a chosen benchmark EoS is very accurate for a wide range of densities and neutron-proton 



asymmetries as well as for large values of the momentum cutoff. At variance with Ref. 12j , we have chosen here the 
SLy5 [l4j Skyrme parametrization to generate the mean-field EoS which has been adopted as a benchmark. 

The article is organized in the following way: in Sec. II the analytical expressions of the second-order contribution 
associated to the full Skyrme interaction are described, and the nature of the divergence is analyzed. In Sec. Ill 
the numerical results for the EoS (with different values of A) are displayed, and the fit of the 9 Skyrme parameters 
is firstly done separately in the case of symmetric (Subsec. III. A), pure neutron (Subsec. III.B), and asymmetric 
(Subsec. IV. C) matter. Finally, we show a global fit for the three cases in Subsec. III.D. Conclusions are drawn in 
Sec. IV. 



II. SECOND-ORDER EQUATION OF STATE WITH THE SKYRME INTERACTION 



As done in Ref. [12|, we treat the EoS of nuclear matter by adding the second-order correction to the first-order 
mean-field energy. The EoS E/A is thus written as the sum of four diagrams which are displayed in Fig. 1 and which 
represent the direct (left) and exchange (right) first-order (upper line) and second-order (lower line) contributions. 
As mentioned in Sec. I, we use in this paper a standard Skyrme interaction like SLy5 in its complete form, namely 

V(n,r 2 ) = to(l + x P„)S(n - f 2 ) + ^(1 + x 1 P a )[k' 2 5(r 1 - r 2 ) + 6(n - r 2 )k 2 } 

+ t 2 (l + x 2 P a )k' ■ 6(n - r 2 )k + \t 3 (l + x 3 P a )p a (R)S(f 1 - f 2 ) 

o 

+ iWofa + <7 2 ) • [k 1 x 6(n - r 2 )k] , (2) 

where R = \(ri + r 2 ), k = ;MVi — V 2 ), k' is the hermitian conjugate of k (acting on the left), P a = 1(1 + ax ■ <j 2 ) 
is the spin-exchange operator. The ten parameters ti, xi, a and Wo characterize a given Skyrme set. In uniform 
matter the spin-orbit force does not play any role, and the associated parameter Wo does not show up in any of the 
quantities discussed below. Consequently, we consider only the remaining nine free parameters. 

Using the general Skyrme force of Eq. ([2]) , the energy per particle in uniform matter can be easily written. The 
mean-field (or HF) result can be found in several papers. The result including both HF and the second-order correction 
reads 
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where 5 is the asymmetry parameter, 



Pn - Pp 
P n + P P ' 
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In Eq. (4) p n and p p are equal to the neutron and proton densities, respectively; 6 = and 5 = 1 represent the 
extreme cases of symmetric and pure neutron matter, respectively. The following notation is used in Eq. (3), 

Gp = \[{l + 8f + {l-5f], 

6 S = 3ti+< 2 (5 + 4a; 2 ), 
O v =t 1 (2 + x 1 ) +t 2 (2 + x 2 ). 

We observe that the second-order term, which is the last term in Eq. j3]), depends not only on p and 8 but on the 
momentum cutoff A as well. Its derivation is discussed in what follows. 

To have a more compact notation, let us write the second-order correction as AE^ by omitting the explicit 
dependence on S, p and A, 

AE^ = (AE + AE x ) nn + (AE + AE X ) PP + AE np . (5) 

The direct AE and the exchange AE X contributions for the neutron-neutron (nn) and proton-proton (pp) channels 
arc written in a box of volume f2 as 
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(6) 
(7) 



where v represents the interaction in momentum space and i denotes either p or n. To simplify the notation, in v 2 (q) 
and v(q)v(k\ — k 2 + q) we have included the factors resulting from the fact that the sum over spin and isospin indices 
has been performed. So these quantities read 

v 2 (q) = 4(7! + 72 ) + (873 + 4 74 )g 2 + 4( 75 + 7 6 )q 4 , 
v(q)v(k 1 -k 2 + q) = (2 7l + 872) + (2 7 3 + 4 74 ) W + (ki - k 2 + q) 2 ] + (875 + 2 76 ) q 2 (k 1 - k 2 + q) 2 . (8) 



The parameters 7$ in Eq. (8) are listed in the following table, 
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to3Xl 2 + X 3h2 


*12 X\2 


*12 + x 12 



where the following notation has been adopted, 

t()3 = t() + TrP", 

6 

, . ^3^3 a 

X03 = t X + — — p a , 
6 

tl2 = tl +t 2 , 
X 12 = tiXi + t 2 X 2 . 

In Eqs. (6) and (7) the energies e are expressed as = , where m* p /m is the effective mass for neutrons or 

protons that we have taken equal to its mean- field value [l4[ , 

771 1 777 1 777 

where u> n = 1 and u p = — 1. In what follows we shall also use the mean- field isoscalar effective mass m*/m, namely 

* \ 1 1 
m„ \ 1 m 



1+ «« e ^' ( 10 ) 
ml 8 Tr 
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and the parameters 



1,1 A 

m* 



mi 



To clarify our compact notation, we stress that in the nn and pp channels, the fci and k 2 momenta appearing in 
the integrals of Eqs. (6) and (7) refer either both to protons or both to neutrons. The integration domain Cj is given 

by C/ = \ki \ < kp;\k 2 \ < kp;\ki + q\ > kp;\k 2 — q\ > kp . This means that k\ and k 2 represent hole states whereas 

k\ + q and k 2 — q represent particle states. The Fermi momentum kp in Cj refers to the proton (neutron) Fermi 
momentum if k\ and k 2 represent both protons (neutrons). 
For the neutron-proton channel one has 



n 3 

(2tt) e 



d 3 fci d 3 k 2 d 3 q 



2v 2 {q) 
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where v 2 {q) is explicitly written in Eq. (8). The factor 2 in Eq. (11) comes from the sum of the np and pn channels. 
The momentum k\ (k 2 ) is associated to a neutron (proton). That means that, in the integration domain Cj that 
we write formally in the same way as for the nn and pp cases, the kp associated to the ki of the neutron (proton) 
represents the neutron (proton) Fermi momentum. We have not specified so far the two different /c_f values in the 
equations to avoid a heavy notation. We will denote later the two Fermi momenta as k n and k p where, 



k„ = ( ^ 2 p n 



3 2 
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We introduce the parameter a depending on k n and k p , 



-n 2 p{l + 5) 



1/3 



1/3 



(12) 
(13) 
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1/3 



< 1. 



(14) 



In general, as already mentioned, one can show that all the corrective terms AE are functions of the density p, 
of the asymmetry parameter 6 and of the momentum cutoff A. Indeed, after some manipulations, the pp and nn 
contributions can be written as the sum of ten terms, 



(15) 



where the first 5 terms describe the nn channel and the last 5 terms describe the pp channel. Let us start with the 
nn case. It can be seen that in the integrals Ij , 



/, (5, p, A) = / du fj (u) [9(1 - u) F( (u) + 0(u - 1) Fi (u)} , j = 1 
Jo 



(16) 



the p, 5 and A dependence enters in the upper limit of integration. The five coefficients Ti and the five functions fi 
are listed in the following table: 



Ti ; 
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r 3 


r 4 


r 5 




/a(«) 


h(u) 


h{u) 


h{u) 


1/15 


1/15 


-i 


4/15 


1/15 


u 


u 2 


u 2 


u A 





5 



The expressions of the ten functions F( (u), F 2 (v) (with j running from 1 to 5) arc provided in Appendix A. The 
coefficients Xj appearing in Eq. (15) (for j = 1 . . . 5) are written as follows 

X2{S,p) = 327r 3 CA£;m*fc r 9 l (to3^i2 + ^03^12), 

Xa{5,P) = 64:Tr 3 C AE m; i kl(t 03 x 1 2 + x 03 t 12 ), (17) 

X5 (5,p) = G^CAEm^it^+x 2 ^), 
where the following notation has been introduced: 

Cae = (18) 

The integral I\ in Eq. (15) has already been encountered in the to — t 3 model treated in Ref. [lH, whereas the 
additional four integrals i2,3.4,5 appear when the full interaction is considered. The last five terms in Eq. (fT5|) (pp 
channel) arc equal to the first five terms with the replacements m* — > m* and k n — > k p in the coefficients x- The 
replacement k n — > k p is also done in the upper limit of the integral in Eq. (16). The coefficients and the functions 
fi and F[ 2 , with i running from 6 to 10, are equal to those already written for i running from 1 to 5. 

For the neutron-proton channel we divide the region of integration into three parts, < |<j| < 2k pi 2k p < \q\ < 2k n 
and \q\ > 2k n and we use the parameter a, defined in Eq. (14). One can derive the following expression, 



AE np IQ^CAEkl m* 



A 15 (bc) s 



r 1 

duuv 2 {2k n u)Ff c (u) + / duuv 2 (2k n u) F 3 abc (u) + / du u v 2 (2k n u) F% bc (u) 

J a Jl 



(19) 

where v 2 (2k n u) = 4(71 + 72) + (873 + 474)(2fc„ u) 2 + 4(75 + 76)(2fc„ u) 4 , and all the other parameters are defined 
above. The expressions of the three functions F? bc (u), Ff c (u) and F$ bc (u) arc provided in Appendix B. 
To summarize, let us write the second-order energy correction in a compact form: 

AE(2 \ X A^ V- AE » + AE u (x , AE n Pfx A x 

— — (6,p,A) = 2^ a.(s,p,A) + —^-(S,p,A) 

i—n,p 
15 

= Y,Xj(S,p)Ij(S,P,A). (20) 
The last 5 terms in the above expression describe the np contributions. For j = 11, . . . , 15 one has 

A 

Ij{5, p, A) = 1 Jj kn du /,-(«) F^{u, p, S), (21) 

where 



fn(u) 




/13(«) 


fu(u) 


/l5(«) 
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u 3 


u 3 


u 5 


u 5 



and 



F^ al {u, P> *) = ^3 [©(« - u)F? bc (u) + Q(u o)0(l - u)Ft c (u) + 0(« - \)Ff c (u)] . (22) 
The last 5 coefficients Xj are equal to 

Xn(S,p) = ^TT 3 C A Em* s kl(tl 3 +xl 3 +t 03 x 03 ), 

Xl2(S,p) = 256TT 3 C A Em* s k n( t 03tl2 + x 03 x 12 ), 

Xi 3 (S,p) = 128TT 3 CA E m:k 9 n (to 3 x 12 +x 03 t 12 ), (23) 

Xu(S,p) = 5127r 3 C , A£;TO*fc ? l l 4i2a;i2, 

Xw(S,p) = $l2<K*CAErKkE{t\ 2 + x\ 2 ). 



6 



Starting from Eq. (20) the asymptotic behavior of the second-order energy contribution can be obtained by taking 
its asymptotic expansion. It can be shown that this leads to 

(6, p, A -»• oo) = 4 p A 5 + 4 p A 3 + al p A + a% p + O ^-f J , (24) 

where the coefficients a 1 depend on 5 and p. We observe that the energy correction diverges as A 5 for large values of 
A and that this divergence (expected from power counting arguments) is much stronger than the linear divergence of 
the to — t 3 model, Eq. (1). 

Finally, the Skyrme EoS up to second order is written as 



E,. 3h 2 /3tt 2 
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A 10m 
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(25) 



III. RESULTS 



A. Symmetric matter (5 — 0) and incompressibility modulus 

The EoS for symmetric matter is given by Eq. (25) by setting (5 = 0. Let us also introduce the pressure and the 
incompressibility modulus calculated up to second order, 



P(p,A) = p 2 ^|(p,A) 



PW(P) 
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Y (x'j V, P)Ij (6, P, A) + X j (8, p)l'j (S, P , A) 



5=0 



(26) 
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Y Xj V, p)Ij (S, p, A) + 2 X ; (S, p% (6, p, A) + xj (5, p)I- (S, p, A) 

i=i 



(27) 



5=0 



where P^ and denote the mean-field (first-order) pressure and the incompressibility modulus, respectively. The 
dependence on the cutoff appears in the second-order corrections. The mean-field expressions depend on the density 
P 0, Ell as follows, 



p (1 Hp) 



lh 2 /3tt 2 \ 2/3 2/3 3 1 /3tt 2 \ 2/ ' 5/3 1 , N a+1 



and 



*S>M = -§£ 

5 m V 2 
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1)*3P' 
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(28) 



(29) 



The ' notation in Eqs. (26) and (27) denotes the derivative with respect to the density p. 

In the upper panel of Fig. 2 we plot the second-order EoS obtained for different values of the cutoff A (see legend) , 
from 0.5 up to 2 fm. -1 . The different equations of state are calculated by using the SLy5 Skyrme parameters and 
are compared with the reference mean-field SLy5 EoS (solid line in (a)). In (b) the second-order correction is plotted 
for the same values of the cutoff. We observe that for a cutoff value equal to 2 fm -1 the correction to the energy at 
the saturation point of nuclear matter, p = 0.16 fm -3 , is very large and amounts to - 80 McV. By using the same 
values of the cutoff A and the SLy5 parameters, the second-order pressure and the second-order incompressibility 
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modulus are displayed in Fig. 3. One may observe how strongly the ultraviolet divergence affects the pressure and 
the incompressibility modulus for large values of the cutoff A. The incompressibility is strongly enhanced by the 
second-order correction and is equal to ~ 625 MeV at the saturation point of matter for A = 2 fm _1 . 

To have a reasonable second-order EoS, we have adjusted the nine parameters of the Skyrme interaction entering 
in the expression of the EoS to reproduce the reference SLy5 mean-field EoS. We have chosen 15 equidistant reference 
points (TV) for densities ranging from 0.02 fm~ 3 to 0.30 fm~ 3 . All the parameters are kept free in the adjustment 
procedure. The minimization has been performed using the following definition for the x 2 



The errors or adopted standard deviations, A_Ej, in Eq. (30) are chosen equal to 1% of the reference SLy5 mean-field 
energies Ei. re f. This choice is arbitrary since we are fitting a theoretical EoS where a standard deviation for this 
quantity has not been estimated. However, the magnitude of the \ 2 defined in Eq. (30) has a clear and reasonable 
meaning: if it is smaller or equal to one, the reference EoS is reproduced within one standard deviation, i.e., within 
a 1% average error by our second-order EoS. The corresponding curves obtained with the adjusted parameters are 
shown in Fig. 4 for different values of A. The quantities which are displayed in this figure are the differences between 
the refitted EoS and the reference SLy5 mean-field EoS for different values of the cutoff A. We observe that the 
deviations are extremely small except at very low densities where they are anyway not larger than 0.06 MeV. In the 
inset of the figure the refitted EoS are plotted and compared with the SLy5 mean-field EoS (solid line). Due to the 
scale, the curves in the inset are practically indistinguishable. The obtained sets of parameters and the x 2 values are 
shown in Table I for each value of the cutoff A. The \ 2 values are always extremely small indicating that, on average, 
the fitted points are deviating much less than 1 % (according to the adopted expression for \ 2 , Eq. (30)) with respect 
to the reference EoS. 

We have noticed that, for the four refitted interactions, the saturation density po and the incompressibility modulus 
are equal in all cases to 0.16 fm~ 3 and 229.9 MeV, respectively. 

The pressure and the incompressibility, Eqs. (|2"6"|) and (|2"T|) . evaluated by using the parameters listed in Table I, 
are plotted in the two panels of Fig. 5. Again, what is plotted is the deviation with respect to the SLy5 mean-field 
reference values. In the two insets, the absolute values are displayed together with the SLy5 mean-field curves (solid 
lines). We stress that the pressure and the incompressibility do not enter in the fits. In spite of this, small deviations 
from the SLy5 reference curves are observed, only at large densities. 



By setting 5 = 1 in Eq. (25) the mean-field plus second-order EoS is obtained for pure neutron matter. The 
ultraviolet divergence with respect to the cutoff is visible in Fig. 6 where the EoS (a) and the second-order correction 
(b) are displayed for different values of the cutoff A. In the upper panel the reference SLy5 mean-field EoS is also 
plotted (solid line). One notices a special and unexpected behavior: starting from the cutoff value A = 1.5 fm -1 
the corrected EoS has an equilibrium point and, for A = 2 fm -1 , the total energy is negative. The appearance of 
an equilibrium point for the second-order EoS of pure neutron matter shows how the ultraviolet divergence is also 
responsible for generating artificial (and unphysical) strong correlations in the system. This anomaly can be cured 
by the adjustment of the parameters. We have performed also in this case the adjustment of the nine parameters of 
the Skyrme interaction with the same definition of \ 2 a s above, Eq. (30). The fitted points are the same as in the 
previous case. In Fig. (7) the deviations with respect to the SLy5 mean-field EoS are shown. Again, the deviations 
are extremely small and are larger at very low densities. In the inset the absolute curves are plotted. The obtained 
parameters and the x 2 values per point are listed in Table II. The % 2 values are extremely small also in this case and 
not larger than 10~ 6 . 



We have chosen the asymmetry value S = 0.5 to illustrate a case of asymmetric matter. The corrected EoS and 
the second-order correction are presented in the upper and lower panels of Fig. 8, respectively. The results of the fit 
(same definition of \ 2 an d same number of fitted points as for the other two cases) are shown in Fig. 9 whereas the 
sets of parameters are listed in Table III. The quality of the fit is very good also in this case as indicated by the low 
values of \ 2 • These values have increased with respect to the two previous adjustments but they still remain much 
lower than 1. 




(30) 



B. Pure neutron matter (S — 1) 



C. An illustration of asymmetric matter (S — 0.5) 
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D. Global fit for the three values 5 = 0, 0.5 and 1 

Finally, a unique and global fit lias been done to readjust the three mean-field plus second-order EoS for symmetric, 
asymmetric and pure neutron matter to reproduce the corresponding SLy5 mean-field curves. The obtained sets of 
parameters are presented in Table IV. In Fig. 10 the three refitted EoS are plotted and in Fig. 11 the deviations from 
the Sly5 mean-field curves are shown. The resulting pressure and incompressibility modulus for symmetric matter 
are shown in Fig. 12. Their deviations with respect to the SLy5 mean-field values are presented in Fig. 13. 

Globally, as one can see from the x 2 values, the quality of the fit is deteriorated with respect to that found for 
each separate case. However, the fit is still of good quality. The \ 2 m this global case is composed by the three 
contributions as calculated in the previous subsections and the final value is divided by three in order to make our 
different results comparable to one another. Specifically, the \ 2 values are still less than 1 up to A = 1 fm _1 . Values 
between 1 and 2 (to be judged by considering the adopted choice of the errors in the expression of x 2 ) are found for 
larger values of the cutoff meaning that the fit is still good. 

We have considered in our global fit only three values of the asymmetry parameter 5, describing symmetric, neutron 
and one case of asymmetric matter, S = 0.5. To judge the quality of the refitted parameters also for other values 
of S (other values of the asymmetry) we show in Figs. 14 and 15 a test performed with the parameters which have 
been obtained from the global fit. Two values of the density (lower and larger than the saturation density) have been 
chosen and the deviations between the refitted EoS and the reference SLy5 mean-field EoS (for different values of the 
cutoff) are plotted as a function of 5 for p = 0.1 (Fig. 14) and 0.2 (Fig. 15) fm~ 3 , respectively. Wc observe that the 
deviations are always reasonably small for all the values of S and (which is the most important result of this test) 
that they do not increase strongly for the values of S which have not been used in the fit. The maximum deviations 
are not larger than 0.4 MeV. 

The values of the saturation density and of the incompressibility modulus for symmetric matter resulting from the 
global fit arc displayed in Table V for the four values of the cutoff A. 

Finally, for the case of the global fits that constitute our more demanding test to the second-order EoS, we have 
estimated the standard deviation of the fitted parameters [l8|]. This analysis allows one to asses how well the used 
reference data toghethcr with the adopted errors constraint the parameters of our model. In particular, the standard 
deviation associated to such parameters are displayed in Table IV. 

IV. CONCLUSIONS 

We have analyzed in this work the nature of the ultraviolet divergence generated by the zero range of the Skyrme 
interaction in the second-order EoS of nuclear matter. The same issue has been previously addressed [l^l but in a 
simple to — ts model and by considering only symmetric nuclear matter. A cutoff regularization has been proposed 
in Ref. [l2| to treat the ultraviolet divergence (which was linear with the momentum cutoff A). In this work the 
velocity-dependent terms of the Skyrme interaction have also been included and both symmetric and asymmetric 
matter have been considered, including the extreme case of pure neutron matter. The expressions of the second-order 
correction to the EoS are derived analytically and a strong divergence (~ A 5 ) is found in the asymptotic expression 
of the corrective terms. A cutoff regularization procedure is adopted first for the single cases of symmetric, neutron 
and asymmetric (5 = 0.5) matter. The resulting fits are of extremely good quality. A global fit is finally performed 
simultaneously for the three EoS. The results are still satisfactory. 

Two interesting conclusions may be drawn: i) Even if the divergence is much stronger than in the simple <o — £3 
case, the fit of the parameters is still possible; ii) The three EoS may be adjusted simultaneously and the problem 
of the appearance of an artificial equilibrium point for neutron matter can be always cured by the adjustment of the 
parameters. 

The adjusted interactions display reasonable properties for nuclear matter. This opens new perspectives for future 
applications of this kind of interactions in beyond mean-field models to treat finite nuclei. It is worth reminding that, 
so far, conventional phcnomenological interactions (adjusted at the mean-field level) have been employed for nuclei in 
different beyond- mean-field calculations [U, [l|| [l7| . 

A drawback of the cutoff regularization procedure is the fact that for each momentum cutoff a different paramctriza- 
tion is generated. A unique set of parameters could be provided by applying the dimensional rcnormalization (men- 
tioned in Sec. I). Work to apply the dimensional renormalization to the second-order EoS of nuclear matter is in 
progress. 
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FIG. 1: Direct and exchange first-order (upper line) and second-order (lower line) contributions to the total energy. 
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FIG. 2: (Color online) (a) Second-order EoS for different values of the cutoff A and (b) second-order correction for symmetric 
nuclear matter calculated with the SLy5 parameters. The SLy5 mean-field EoS is also plotted in (a) (solid line). 



TABLE I: Parameter sets obtained in the fit of the EoS of symmetric matter for different values of the cutoff A compared with 
the original set SLy5. In the last column the \ 2 values are shown. 



xi 



■r-2 



•I'M 



to ti £2 £3 a;o 

(MeV fm 3 ) (MeV fm 5 ) (MeV fm 5 ) (MeV fm 3+3a ) 
SLy5 -2484.88 483.13 -549.40 13763.0 0.778 -0.328 -1.0 1.267 0.16667 



A(fm" 1 ) 

0.5 -1817.280 

1.0 -1132.001 

1.5 -608.125 

2.0 -331.658 



646.948 4373.135 

807.361 -323.413 

71.647 241.517 

660.677 -695.979 



10101.307 -0.0002 -3.464 -1.314 6.233 

7555.400 0.733 1.201 0.644 5.012 

-3920.616 1.565 -2.376 1.655 6.111 

-90.060 3.000 -0.803 -1.120 164.031 0.754 3.3e-06 



X 

0.246 6.2e-06 
0.457 4.9e-07 
0.834 7.2e-06 
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P (fm 3 ) 

FIG. 3: (Color online) Second-order pressure (a) and incompressibility modulus (b) calculated with the SLy5 parameters for 
different values of the cutoff. The mean-field SLy5 curves are also plotted in both panels (solid lines). 

TABLE II: Parameter sets obtained in the fit of the EoS of neutron matter for different values of the cutoff A compared with 
the original set SLy5. In the last column the \ 2 values are shown. 



to tl t% t'i Xo Xl X2 xz a 

(MeVfm 3 ) (MeVfm 5 ) (MeV fm 5 ) (MeV fm 3+3Q ) 

SLy5 -2484.88 483.13 -549.40 13736.0 0.778 -0.328 -1.0 1.267 0.16667 

A(fm- 1 ) x 2 

0.5 -535.222 403.303 1660.746 42905.115 0.094 -0.970 -1.031 1.094 0.144 9.1e-09 

1.0 -1941.276 92.989 393.422 -137583.116 0.609 -0.502 -1.010 1.057 0.613 1.5e-06 

1.5 -18033.283 319.198 -186.907 110184.232 1.846 -1.113 -0.929 1.893 0.006 4.7e-08 

2.0 -218.464 598.755 -538.604 496.206 0.015 -0.885 -0.745 14.793 0.205 6.6e-06 



TABLE III: Parameter sets obtained in the fit of the EoS of asymmetric matter (S = 0.5) matter for different values of the 
cutoff A compared with the original set SLy5. In the last column the \ 2 values are shown. 
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-0.328 


-1.0 
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0.16667 




A'fm" 1 ) 
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-275.173 
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-0.434 
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FIG. 4: (Color online) Deviations between the refitted second-order EoS (for different values of the cutoff) and the SLy5 mean- 
field curve for symmetric nuclear matter. In the inset the absolute values are plotted and compared with the SLy5 mean-field 
EoS (solid line). 



TABLE IV: Parameter sets obtained in the fit of the EoS of symmetric, asymmetric and pure neutron matter for different 
values of the cutoff A compared with the original set SLy5. The standard deviation, a, estimated for the different parameters 
is also given. In the last column the % 2 values are shown. 





to 


h 


ti 


ts 


Xo 


Xl 


X2 


x 3 


a 






<Tt 




(Tt 2 


a t3 


(Txo 














(MeV fm 3 ) 


(MeV fm 5 ) 
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SLy5 


-2484.88 


483.13 


-549.40 


13736.0 


0.778 


-0.328 


-1.0 


1.267 


0.16667 




A(fm" 1 ) 




















x 2 


0.5 


-2022.142 


290.312 
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12334.459 


0.481 


-5.390 


-1.304 


0.880 


0.259 


0.411 
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0.01415 
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0.045037 


0.000784 




2.0 


-718.397 


573.884 


-497.766 


6179.243 


0.391 


-0.393 


-0.574 


0.785 


1.051 


1.313 




0.343 


0.251 


0.261 


8.33 
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0.001850 
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FIG. 6: (Color online) As in Fig. 2 but for pure neutron matter. 
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FIG. 8: (Color online) As in Fig. 2 but for asymmetric nuclear matter in the case 8 = 0.5. 
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TABLE V: Saturation density and incompressibility modulus resulting from the global fit for symmetric nuclear matter. 



A (fm" 1 ) 


po (fm" 3 ) 


Koo (MeV) 


0.5 


0.16 


236.36 


1.0 


0.16 


230.52 


1.5 


0.16 


236.28 


2.0 


0.16 


222.76 
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V. APPENDIX A 



The expressions of the ten functions F((u) and F^iu) (with j running from 1 to 5) appearing in Eq. (16) read, 



Fl(u) 
F 2 (u) 

F?(u) 



F?(u) 



15 



15 



4 + —u - 5m 3 + -m 5 log(l + u) + 4 - — u + 5m 3 - -m 5 log(l - u) + 29m 2 - 3m 4 - 40m 2 log 2 



(4 - 20m 2 - 20m 3 + 4m 5 ) log(u + 1) + (20m 2 - 20m 3 + 4m 5 - 4) log(u - 1) + 22m + 4m 3 + (40m 3 - 8m 5 ) log u. 



15 _ _6_ 
~8 ~ 7u 



12m- 



195 2 129 4 267 6 \ , 

u 2 + — u 4 -—u 6 ) log (I- u) 



311m 1457m 3 



267m 5 



28 



14 



2S 



/15 6 „ 195 2 129 4 267 fi \, „ , , _ , „ 

- f -g + Tjr - 12« - — + — - -^w 6 J log(l + u) + (8m - 140u 3 ) log 2 



58 550u 2 



92m 4 



6 92 
— - 16m + 70m 3 - 68m 4 + yM 6 ) log(M - 1) 



184 



7 



6 



92 



136m 4 - — u° logM + -— + 16m - 70m 3 - 68m 4 + —u° \og{u + 1) 



7u 



15 12 „ 135 2 89 4 183 A, , 

h 16m m 2 H m 4 m log(l - m) 

4 7m 4 4 28 J BV ' 



7 



311m 1051m 3 183m 5 



14 



7 



14 



H + B - 16m - H^m 2 + ^M 4 - ^m 6 ) log(l + u) + (16m - 200m 3 ) log 2 



Fi(u) = — 



116 792m 2 128m 4 

tt + + 
7 7 7 



12 128 

24m + 100m 3 - 96m 4 + m 6 I log(u - 1) 

7m 7 



192m 4 



15 12 



256 



7 



12 



128 



u° logM + h 24m - 100m 3 - 96m 4 H M b log(u + 1) 



7m 



45 



29 



57 



1- 4m ir H m 4 M b log(l - it) - 



7u 



28 



7 



311m 442m 3 57m 5 



14 



7 



14 



(l6M-80M 3 )log2 



15 12 , 45 2 29 4 57 6 . 
-- - + 4m + yM 2 - yM 4 + ^ ) log(l + u) 



F%{u) = — 



116 330m 2 44m 4 

tt + + 
7 7 7 



— - 12m + 40m 3 - 36m 4 + — m 6 I log(u - 1) 
7m 7 



72m 4 - yM b ) logM+ ( -— + 12m-40m j 



311m 239m j 



420 



210 
2 



it 
28 



12 

7m 

8 4 5 



1 



8 35^ ~ 8 U 
8m 4 
105 



40 



56 



log2-[i 
log[l + m] 



36m 4 

2 L 2 
35u ~ 8 U 



44 

— M 6 I bg(M + 1) 



40 



56 



log[l - u] 



_ B/ , 58 88m 2 

F 9 5 (m) = 

2y ' 105 105 105 



"16 4 


16 - 
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logM — 


M - 


15 


105 





2 4 2 o 8 4 

M + -M M 

35m 15 3 15 



105 



log[M — 1] 



35m 15 



15 



105 
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VI. APPENDIX B 

The expressions of the functions Fi bc (u), Ff bc (it) and Fg bc (u) appearing in Eq. (19) are 
ii 2 

F? bc (u) = — (ll6 4 c + I8ab 3 c 2 + I8a 2 b 2 c 3 + lla 3 6c 4 ) + w 3 (186 3 c - 96 4 c + 6a6 2 c 2 - 96 3 c 2 - 3a6 3 c 2 + 18a 2 6c 3 
-3a6 2 c 3 - 9a 2 6 2 c 3 - 9a 2 6c 4 ) + u 4 (66 2 c - 66 3 c + 26 4 c + 6a6c 2 - 66 2 c 2 - 6ab 2 c 2 + 36 3 c 2 - \ab 3 c 2 
-6abc 3 - h 2 c 3 + 3ab 2 c 3 + 2a6c 4 ) + u 5 (4bc - 66 2 c - b 3 c + -6 4 c - 66c 2 + 66 2 c 2 + h 3 c 2 - be 3 

+ \b 2 c 3 + ^6c 4 ) 

+ log 2 1 - 20u 2 (6 4 c + a 3 6c 4 ) + u 3 (306 4 c - 606 3 c + 306 3 c 2 - 60a 2 6c 3 + 30a 2 6 2 c 3 

+30a 2 6c 4 ) + u 5 (806c - 1206 2 c + 806 3 c - 206 4 c - 1206c 2 + 1206 2 c 2 - 406 3 c 2 + 806c 3 - 406 2 c 3 

-206c 4 ) | + 10u 5 6c(2 - 6 - c) 3 log(2(2 - 6 - c)u) 

+ | i(6 5 - 5aW - 5a 3 6 2 c 3 + a 5 c 5 ) + I^(6 4 - 6 5 - 2a W + 2a 2 6 3 c 2 + a 4 c 4 - a 4 6c 4 ) 
+10u 2 (6 3 - 26 4 + 6 5 + a 3 c 3 - 2a 3 6c 3 + a 3 6 2 c 3 ) + 10u 3 (6 2 - 36 3 + 36 4 - 6 5 + a 2 c 2 - 3a 2 bc 2 + 3a 2 6 2 
-a 2 6 3 c 2 ) + it 5 (-4 + 206 - 406 2 + 406 3 - 206 4 + 46 5 )| log (6 + ac + 2u - 2bu) 

+ |(i^ a 4 &c 4 u + i 0a V(26c 3 - 6 2 c 3 ) + aV(306c 2 - 306 2 c 2 + y6 3 c 2 ) 

+u 5 (-206 + 406 2 - 306 3 + 106 4 - ^6 5 ) j log(2ac + 4u - 2bu) 

+ | i(6 5 - 5a W - 5aW + a 5 c 5 ) + ^(6 4 - 6 4 c - 2aW + 2aW + a 4 c 4 - a 4 c 5 ) 
+10u 2 (6 3 - 26 3 c + 6 3 c 2 + a 3 c 3 - 2a 3 c 4 + a 3 c 5 ) + 10u 3 (6 2 - 36 2 c + a 2 c 2 + 3b 2 c 2 - 3a 2 c 3 - 6 2 c 3 
+3a 2 c 4 - a 2 c 5 ) + w 5 (-4 + 20c - 40c 2 + 40c 3 - 20c 4 + 4c 5 )| log(6 + ac + 2u - 2cw) 

{15 15 
— 6 4 cu + w 2 (206 3 c - 106 3 c 2 ) + u 3 (306 2 c - 306 2 c 2 + y& 2 c 3 ) + u 5 (40c 2 - 20c - 30c 3 + 10c 4 

-^c 5 )|log(26 + 4u-2cu) 

J I(_ 6 6 + 5a 2 6 3 c 2 + 5a 3 fe 2 c 3 _ ^5) + i^(_ fe 4 + & 5 + ^ + ^2 Q 2 _ ^2^2 _ ^2^3 

1 2 4 

-a 4 c 4 + a 4 6c 4 + a 4 c 5 ) + u 2 (-106 3 + 206 4 - 106 5 + 206 3 c - 206 4 c - 106 3 c 2 - 10a 3 c 3 + 20a 3 6c 3 
-10a 3 6 2 c 3 + 20a 3 c 4 - 20a 3 6c 4 - 10a 3 c 5 ) 

+it 3 (-106 2 + 306 3 - 306 4 + 106 5 + 306 2 c - 606 3 c + 306 4 c - 10a 2 c 2 + 30a 2 6c 2 - 306 2 c 2 - 30a 2 6 2 c 2 
+306 3 c 2 + 10a 2 6 3 c 2 + 30a 2 c 3 - 60a 2 6c 3 + 106 2 c 3 + 30a 2 6 2 c 3 - 30a 2 c 4 + 30a 2 6c 4 + 10a 2 c 5 ) 
+it 5 (4 - 206 + 406 2 - 406 3 + 206 4 - 46 5 - 20c + 806c - 1206 2 c + 806 3 c - 206 4 c + 40c 2 - 1206c 2 

+1206 2 c 2 - 406 3 c 2 - 40c 3 + 806c 3 - 406 2 c 3 + 20c 4 - 206c 4 - 4c 5 ) I log (6 + ac + 2u - 2bu - 2cu) 
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+ { - ^-a 4 bc 4 u + u 2 (-20a 3 bc 3 + I0a 3 b 2 c 3 + 20a 3 bc 4 ) + u 3 (-30a 2 bc 2 + 30a 2 6 2 c 2 - ^-a 2 b 3 c 2 + 60a 2 bc 3 
4 2 

-30a 2 6 2 c 3 - 30a 2 6c 4 ) + u 5 (206 - 40b 2 + 30b 3 - 10b 4 + ^6 5 - 80bc + 1206 2 c - 606 3 c + 106 4 c + 1206c 2 
-1206 2 c 2 + 306 3 c 2 - 806c 3 + 406 2 c 3 + 206c 4 ) j log(2ac + Au - 2bu - Acu) 

+ j - ^6 4 cu + u 2 (-206 3 c + 206 4 c + 106 3 c 2 ) + u 3 (-306 2 c + 606 3 c - 306 4 c + 306 2 c 2 - 306 3 c 2 - ^6 2 c 3 ) 
+u 5 (20c - 806c + 1206 2 c - 806 3 c + 206 4 c - 40c 2 + 1206c 2 - 1206 2 c 2 + 406 3 c 2 + 30c 3 - 606c 3 + 306 2 c 3 - 10c 
+ 106c 4 + ^c 5 ) j log(4u - 2cm + 26 - 46u) 

_ J 1( 6 5 _ 5a 2 6 3 c 2 _ 5a 3 6 2 c 3 + ^5) + _ ^2^2 + ^ + 10u 2 (fe 3 + ^3) 



+ 10u 3 (6 2 + a 2 c 2 ) - 4m 5 j log(6 + ac + 2u) 

Ff c {u) = {^(b 5 -5a 2 b 3 c 2 -5a 3 b 2 c 3 + a 5 c 5 ) + ^(-b 4 + 2a 2 b 2 c 2 -a 4 c 4 ) 
+10w 2 (6 3 + a 3 c 3 ) - 10w 3 (6 2 + a 2 c 2 ) + 4u 5 } log(2u - ac - 6) 

+ {I( fo 5 _ 5a 2 5 3 c 2 + 5a 3 6 2 c 3 _ ^5) + l^ (fe 4 _ ^2^2 + ^4) 

2 4 
+ 10u 2 (6 3 - a 3 c 3 ) + 10u 3 (6 2 + a 2 c 2 ) - 4u 5 } log(2u - ac + 6) 

+{~(-b 5 + 5a 2 b 3 c 2 - 5a 3 b 2 c 3 + a 5 c 5 ) + — (6 4 - 2a 2 6 2 c 2 + a 4 c 4 ) 
L 2 4 

+10w 2 (-6 3 + a 3 c 3 ) + 10w 3 (6 2 + a 2 c 2 ) - Au 5 } log(2u + ac - 6) 

-{V - 5aW - 5a 3 6 2 c 3 + a 5 c 5 ) + — (6 4 - 2a 2 6 2 c 2 + a 4 c 4 ) 
2 4 

+ 10u 2 (6 3 + a 3 c 3 ) + 10w 3 (6 2 + a 2 c 2 ) - 4u 5 } log(2u + ac + 6) 
+ lla6cu(6 2 + a 2 c 2 ) + 4a6cu 3 
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Ff c (u) = — (ab 4 c + a 3 b 2 c 3 ) + u 2 [18ab 3 c - 9a6 4 c + lla 3 6c 3 - — a 3 6 2 c 3 ) + u 3 (6ab 2 c - 6ab 3 c + 2ab A c) 

+ {-h 5 + \a 2 b 3 c 2 - \a 3 b 2 c 3 + ia 5 c 5 + u(-^-b 4 + + ^aW - ^aW - ^a 4 c 4 

+ ^a 4 6c 4 ) + u 2 {-10b 3 + 20b 4 - 106 5 + 10a 3 c 3 - 20a 3 oc 3 + 10a 3 6 2 c 3 ) 
+?i 3 (-106 2 + 30o 3 - 30o 4 + 10o 5 - 10a 2 c 2 + 30a 2 oc 2 - 30a 2 o 2 c 2 + 10a 2 b 3 c 2 ) 

+w 5 (4 - 206 + 40b 2 - 406 3 + 20o 4 - 4b 5 ) | log(6 - ac + 2u - 2ou) 

+ [ IftS _ 5 2 5 3 c 2 _ 5 3 6 2 c 3 + 1 + _ ^5 _ ^2^2 + 2^ + _ ^4^4 

I _■ TC Tt jfi _ 'rt ~ t 

+u 2 (106 3 - 20o 4 + 10b 5 + 10a 3 c 3 - 20a 3 6c 3 + 10a 3 6 2 c 3 ) 

+w 3 (10o 2 - 30b 3 + 306 4 - 106 5 + 10a 2 c 2 - 30a 2 oc 2 + 30aVc 2 - 10a 2 o 3 c 2 ) 

+u 5 (-4 + 206 - 40b 2 + 40b 3 - 206 4 + 4o 5 j log(6 + ac + 2u - 2bu) 
+u 5 (~20b + 40b 2 - 30b 3 + 106 4 - ^6 5 ) j log(2ac + 4u - 2bu) 

J h l _ ^ a 2 fe 3 c 2 _ 5 a 3 & 2 c 3 + £^ +u {— - — aW + — a 4 c 4 ) +u 2 (106 3 + 10a 3 c 3 ) 
122 2 2 42 4 

+ M 3 (106 2 + 10a 2 c 2 ) - 4u 5 | log(6 + ac + 2u) 

+-(& - ac + 2u) 3 (46 2 + 4a 2 c 2 - 6acu - 4m 2 + 36(4ac + 2m)) log (6 - ac + 2u) 
8 

+ — [-6ac - 4u + 2ou] [2ac - 4u + 2ou] 3 log(4w - 2ow - 2ac) 



